Double beta decay versus cosmology: Majorana CP 
phases and nuclear matrix elements 



Frank Deppisch 1 *, Heinrich Pas 1 ^ Jouni Suhonen 2 * 

1 Institut fur Theoretische Physik und Astrophysik, 
Universitdt Wiirzburg, D-97074 Wiirzburg, Germany 

2 Department of Physics, University of Jyvdskyld, 

P.O.B. 35, FIN-40014, Jyvdskyld, Finland 

Abstract 

We discuss the relation between the absolute neutrino mass scale, the effective mass 
measured in neutrinoless double beta decay, and the Majorana CP phases. Emphasis 
is placed on estimating the upper bound on the nuclear matrix element entering 
calculations of the double beta decay half life. Consequently, one of the Majorana 
CP phases can be constrained when combining the claimed evidence for neutrinoless 
double beta decay with the neutrino mass bound from cosmology. 
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Over the past years much effort has been invested to probe leptonic mixing with increasing 
accuracy, and, in fact, a unique picture is evolving from the precise measurements of neutrino 
oscillation probabilities. For a full construction of the mixing matrix, however, knowledge 
about CP violating phases is necessary. 

To describe leptonic mixing one can always work in a basis where the charged lepton Yukawa 
matrix is diagonal. In this case the neutrino mass matrix m u can be written in the flavor 
basis as 

= U*m di * s U\ (1) 

with m diag = diag(mi, m 2 , m 3 ), and the Maki-Nakagawa-Sakata (MNS) matrix U, which 
contains three mixing angles and one CP violating Dirac phase. The determination of the 
mixing angles is subject to neutrino oscillation experiments, as is (at least in principle) the 
determination of the Dirac phase - the leptonic analogue of the CKM phase in the quark 
sector. Three of the mixing angles can be identified with the maximal, large and small 
observables measured in atmospheric, solar and reactor neutrino oscillations, respectively 
[1]. Important information on the Dirac phase can be expected from a combination of future 
long-baseline experiments [2], and ultimately, from a neutrino factory [3]. 
If the neutrino is of Majorana type, two more phases enter, though. The determination 
of these Majorana phases is, in fact, the most challenging task in the reconstruction of 
the fundamental parameters of the Standard Model particle content. In general, U can be 
written as 

U = V ■ diag(l, e^ 12/2 , e i<te/2 ), (2) 
where V is parametrized in the standard CKM form, 



V 
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(3) 



and (fiij are the Majorana phases under discussion. A viable possibility to obtain information 
on the Majorana phases is to compare measurements of the absolute neutrino masses 
with the elements of the neutrino mass matrix m" in the flavor basis. While absolute 
neutrino masses are most stringently constrained from cosmology, only the ee element of m v 
is experimentally accessible, being the effective mass m ee measured in neutrinoless double 
beta decay. Several works have discussed the relations of m ee , absolute neutrino masses 
and Majorana phases [4]. In [5] it was pointed out that one could restrict the Majorana 
Phase 12 by using the recently claimed evidence for neutrinoless double beta decay [6] and 
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the cosmological neutrino mass bound derived by the WMAP collaboration [7] , if a certain 
nuclear matrix element (NME) calculation [8] is assumed. An important issue is, however, 
the uncertainty in the neutrino mass determination within the double beta decay framework 
due to systematical limitations in such NME calculations. In this work we focus on what 
can be learned about Majorana phases from recent double beta decay and cosmological 
structure formation data, in view of an upper bound on the NME. In the following, we 
first review the claimed evidence for neutrinoless double beta decay and discuss the upper 
bound on the NME. Finally, we compare the lower bounds on m ee obtained with the upper 
bounds on neutrino masses from cosmology. 
The half life of neutrinoless double beta decay is given by 
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where m e denotes the electron rest mass, is a phase space factor, and the NME is 

given by M {0u) =M G t-M f , being a combination of Gamov- Teller and Fermi transitions. 
The neutrinoless double beta decay is sensitive to the ee element of the mass matrix m v (Q) 
in flavor space, 



m. 
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23- 



Here the mass eigenstates are expressed as mi, m 2 = \Jm 2 + Am 2 2 , m 3 = 
where the plus (minus) sign applies in the normal (inverted) mass hierarchy case. 
From the Chooz and Palo Verde experiments we know that \V^ 3 \ <C \V^\, \ V^ 2 \- Moreover, 
in the quasi-degenerate mass range explored by the present experiments, one has 3> 
Amf 2 , Am 23 . Employing the above approximations we arrive at a very simplified expression 
for m ee (see, e.g. [9]): 



1 - sin 2 (2# 12 ) sin 2 
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illustrating that m ee is mostly sensitive to #12 and 0i2- 

In Figs, ^and m\/\m ee \ 2 is shown as a function of 0i2, using the exact relation (jSJ). The 
full colored/shaded region indicates the allowed range according to the combination of 2a 
(Fig. |U) and 3a (Fig. 01) limits on the neutrino oscillation observables #12, 6*23, #13, Am\ 2l 
Am 23 [1]. The remaining Majorana phase 023 is varied in its full range, 023 £ [0,27r]. The 
dark/red bands correspond to a fixed maximal (upper band), best-fit (middle) and minimal 
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Figure 1: mf/\m ee \ 2 as a function of the Majorana phase 0i2- The full colored region is 
the allowed range defined by present 2a limits on 6 I 12 , 2 3, 6> 13 Amf 2 , Am^. The dark/red 
bands correspond to a fixed maximal (upper band), best-fit (middle) and minimal (lower) 
value for 9u, varying all other parameters. The second Majorana phase $23 has been varied 
in the full range [0, 2%}. 

(lower) value for # 12 , varying all other parameters. The minimum values for m\/\m ee \ 2 at 
012 = 7r are at 3.81 and 3.23 for 2a and 3a oscillation limits, respectively. Thus, if the 
experimental values for m\/\m ee \ 2 turn out to be smaller, a Majorana CP phase 0i 2 = vr is 
excluded. If, finally m\/\m ee \ 2 < 1, the bound from unitarity of the MNS matrix U would 
be violated, resulting in the conclusion, that either the limit on mi or the limit on \m ee \ is 
not applicable. 

Recently, a new publication of data of the Heidelberg- Moscow experiment [10], searching 
for the double beta decay of 76 Ge, has appeared. In this work the authors have analyzed 
the data taken in the period 1990-2003, and applied a new energy calibration, which in- 
creased the previous claim for evidence [6] to 4.2a statistical significance. The allowed range 
corresponds to [11] 

\m ee \ =f.(0.39 - 0.49 eV) [la], (7) 
\m ee \ = f. (0.32 - 0.54 eV) [2a], (8) 
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Figure 2: As Fig. 1, but with 3a limits on 6 12 , #23, 6*i 3 Am^ 2 , Am^. 

\m ee \ = £• (0.24 - 0.58 eV) [3a]. (9) 

Here f = M SMK /M {0u) denotes the normalization to the NME M SMK = 4.2, calculated 
in the pn-QRPA model of [8]. While the initial claim caused a critical debate [12] and 
was not confirmed by an independent analysis of the data [13], several of these issues have 
been clarified in [14]. Remaining criticisms concern the exact peak position in the energy 
spectrum and the relative strength of measured background lines, which could suggest, 
that the observed signal is due to an unidentified background line or that the statistical 
significance of the signal is overestimated [15]. In any case, we feel motivated to take the 
evidence claim at face value and discuss possible consequences, although we stress that an 
independent test of the claimed evidence is essential, if possible with a different double beta 
emitter isotope. Such a test could be realized by the recently started CUORICINO [16] 
and NEMO experiments [17] and the recent MPI proposal [18] which revived the GENIUS 
proposal of the Heidelberg group [19]. 

Any conclusions about relations of the absolute neutrino mass mi and the double beta 
decay observable m ee depend crucially on the magnitude of the calculated NME. Typically 
the uncertainty of such calculations has been estimated to be a factor 2-3, assumed around 
a given central value such as the calculation in [8], i.e. £ G [0.5,2]. In the following we 
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will argue in favor of a more stringent upper bound on £ and that it constrains the allowed 
range of the CP Majorana phase 4> 12 . In Table 1 a scan of all the available matrix element 
calculations has been performed. Most of the used models are based on the proton-neutron 
quasiparticle random-phase approximation (pnQRPA), like the renormalized pnQRPA, de- 
noted by RQRPA in the table, the self-consistent pnQRPA (SQRPA), the self-consistent 
RQRPA (SRQRPA), and the fully self-consistent RQRPA (full-RQRPA) . In addition, the 
pnQRPA has been improved by performing a particle- number projection on it in Ref. [20]. 
This theory has been denoted by projected pnQRPA in the table. The line with pnQRPA 
+ pn pairing in the table denotes a theory where proton-neutron pairing has been added 
to the RQRPA framework. 

In these calculations various sizes of the proton and neutron single-particle valence spaces 
have been used. They range from rather modest to very extensive single-particle bases. 
Also, the single-particle energies have been obtained either from the experimental data, 
or, more frequently, from a Coulomb-corrected phenomenological Woods-Saxon potential 
where the parameters have been adjusted to reproduce spectroscopic properties of nuclei 
close to the beta-stability line. In addition, in some calculations the Woods-Saxon single- 
particle energies have been varied close to the proton and neutron Fermi levels to reproduce 
low-energy spectra of the neighboring nuclei with odd number of protons or neutrons. Hence, 
remarkably diversified starting points have been used for the calculations. 
Some nuclear matrix elements of this table can be discarded due to various deficiencies in the 
theoretical frameworks used to evaluate them. This concerns the shell-model matrix element 
j^(0u) _ 5 go f Haxton and Stephenson [21], who used the weak-coupling approximation 
in evaluation of it. This is quite a rough approximation, and the more recent matrix 
element = 1.74, obtained by performing a large-scale shell- model calculation with 

realistic two-body forces, should be more reliable, although some doubts concerning the 
adequacy of the size of the used single-particle basis have been voiced. The reliability of 
the pnQRPA calculation including the proton-neutron pairing has been questioned due to 
the way the pairing is introduced to the theory. One can ignore the corresponding matrix 
element if one wants, without changing our final conclusions concerning the upper limit of 
the computed NME. The largest matrix element was calculated by Tomoda et al. in [22] 
by using a quasiparticle mean-field based VAMPIR approach with particle-number and 
angular-momentum projections included. The shortcoming of this approach is that it does 
not include the proton-neutron interaction in its framework, being essential in description 
of charge-changing nuclear transitions, which occur also in the double beta decay. 
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The above considerations lead to the range 

0.59 > M (0iy) > 4.59 (10) 

for the acceptable nuclear matrix elements. Since the upper limit of the above range comes 
from a pnQRPA calculation, one may ask how does the u g pp problem" of the pnQRPA 
affect this value. This problem concerns the calculated matrix element of the two-neutrino 
double beta decay Ai^, which turns out to depend strongly on the parameter g pp , used as 
a scaling parameter of the particle-particle part of the proton-neutron two-body interaction 
[23,24]. However, while the uncertainty in g pp affects the lower bound on NME calculations 
dramatically (Ai^ can even become zero for a large value of g pp ) towards lower values of 
g pp the value for the NME enters a plateau, making the upper bound more stable against 
variations in g pp . Moreover, even though the NME corresponding to the two-neutrino 
mode depends strongly on g pp within its physical range, the NME corresponding to the 
neutrinoless mode depends only very weakly on this parameter. This can be clearly seen in 
Figs. 1 and 2 of Ref . [20] , where the relevant double Gamow-Teller and double Fermi matrix 
elements have been plotted as functions of g pp for the pnQRPA and projected pnQRPA 
calculations of the 76 Ge double beta decay. The variation of these matrix elements around 
the physical value of g pp ~ 1 is less than 20 per cent. Adding an uncertainty of 20 per cent 
to the above range of acceptable values of NMEs leads us to the upper limit 

M {0u) < 5.5 (11) 

of the NME. 

Consequently, the limits (0)-© read as 

|m ee |>0.30eV [la], (12) 
|m ee |>0.24eV [2a], (13) 
|m ee | > 0.18 eV [3a]. (14) 

These lower bounds on m ee have to be compared to the most stringent upper bounds on 
the absolute neutrino mass scale, mi, which are presently provided by data on cosmological 
structure formation. According to Big Bang cosmology, the masses of nonrelativistic neutri- 
nos are related to the neutrino fraction of the closure density by J2i m i — 40 Q u /ig 5 eV, where 
h e5 is the present Hubble parameter in units of 65 km/(s Mpc). In the currently favored 
ACDM cosmology with a non-vanishing cosmological constant A, there is scant room left 
for the neutrino component. The free-streaming relativistic neutrinos suppress the growth 
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NME 


Theory 


References 


5.00,1.74 


Shell model 


[21,25] 


1.53-4.59 


pnQRPA 


[8,20,26-33] 


1.50 


pnQRPA+ pn pairing 


[28] 


3.45 


projected pnQRPA 


[20] 


6.76 


VAMPIR 


[22] 


1.87-2.81 


RQRPA 


[32-35] 


0.59-0.65 


SRQRPA 


[32] 


2.40 


full-RQRPA 


[33] 


3.21 


SQRPA 


[33] 



Table 1: Compilation of calculated nuclear matrix elements for the neutrinoless double 
beta decay of 76 Ge. The first column gives the value(s) of the NME, the second column the 
theory used to evaluate the NME, and the last column the works where the quoted theory 
was used to evaluate the NME. 

of fluctuations on scales below the horizon (approximately the Hubble size c/H(z)) until 
they become nonrelativistic at redshifts z ~ rrij/3T ~ 1000 (rrij/eV). 
Recent limits, obtained by combining CMB measurements with data on the large scale 
structure of the universe, imply an upper bound on the sum of the three neutrino mass 
eigenstates [36], 

J2m < 0.42 - 1.80 eV [2a]. (15) 

i 

The exact cosmological bound depends on the data and specific priors used in the analysis: 
The weakest bound utilizes only data from WMAP and the Sloan Digital Sky Survey. In 
particular, limits utilizing the Lyman-o; forest, the absorption observed in quasar spectra 
by neutral hydrogen in the intergalactic medium, provide more stringent bounds [37] 

m < 0.42 - 0.69 eV [2a], (16) 

i 

as compared to the data sets without the Lyman-a forest [38], 

Y^rrii < 0.6 - 1.8 eV [2a]. (17) 

i 

Since the effects of possible systematics in this data set need still to be explored further, 
both observationally and theoretically, each set of analyses should be discussed here. Within 
the present decade, the combination of SDSS data with CMB data of the PLANCK satellite 
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Figure 3: Relation of the cosmologically relevant sum of neutrino masses, J2i m i an d the 
lightest neutrino mass mi. The upper and lower branches correspond to normal and inverse 
hierarchy, respectively. The light/blue and dark/red bands indicate the 2a and 3a range of 
Am| 2 and Am^, while the central line marks the best fit. 

will obtain a 2 a detection threshold on J2i m i close to 0.1-0.2 eV [36,39]. Fig. |3]shows the 
relation of the sum of neutrino masses with mi, 



vrii = m 1 1 + 



\ 



Am 2 
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mf 



+ 



\ 



1 Amj 2 ± Ami 3 
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rrii 



2 

mf 



;is) 



for normal and inverse hierarchy in the upper and lower panel, respectively. The different 
curves indicate the best-fit and upper and lower 2a and 3a ranges. 

Combining the resulting bound on m x with the range for the effective mass m ee given in [10], 
one obtains: 

m\/\m ee \ 2 < 0.36 - 0.93 (0.64 - 1.7) at 2a (3a) (19) 
for data sets using the Ly-a forest and 



m\/\m ee \ 2 < 0.71 - 6.3 (1.3- 11.1) at 2a (3a) 



(20) 
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without the Ly-a forest. It is obvious, that the 2a range for the double beta decay observable 
is in conflict with all cosmological fits including the Ly-a forest and with some without 
the Ly-a forest. This may indicate that either the double beta decay signal is due to a 
statistical fluctuation or due to a mechanism involving exchange of other particles besides 
light massive Majorana neutrinos. Examples for the latter case include, e.g. sparticles in 
i?-parity violating supersymmetry, leptoquarks, or right-handed neutrinos and W bosons 
(for an overview see [40]). Alternatively, the cosmological neutrino mass bound may not be 
applicable, e.g. by introducing broken scale-invariance in the primordial power spectrum 
[41], or due to fast decays of the relic neutrino background [42]. The 3a range for m ee , 
however, is still compatible with most of the cosmological bounds although in most cases 
it is smaller than the minimum value, m\/\m ee \ 2 = 3.81 for <pi2 — tt. If a value 1 < 
m\/\m ee \ 2 < 3.81 will be confirmed in future cosmological data fits or in the upcoming 
tritium beta decay spectrometer KATRIN [43], Majorana CP phases around 0i2 = vr can 
be excluded. 

In conclusion, we discussed upper bounds on nuclear matrix elements and their implications 
for the Majorana CP phase 12 , when the recent evidence claim for neutrinoless double beta 
decay and the cosmological neutrino mass bound are combined. We deduced that for a 
combination of 2a experimental limits in most analyses the mass mechanism interpretation 
of the double beta decay evidence is incompatible with the cosmological neutrino mass 
bound. On the other hand, the range of the Majorana phase </>i 2 can be constrained, when 
combining 3a experimental limits with reasonable upper bounds for the nuclear matrix 
elements. Assuming CP conservation, the CP phase factors exp(i0jj) are reduced to CP 
parities rjij = ±1. In this case the neutrino CP parity is fixed to 7712 = +1. 
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